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Abstract 

J> I We review the current status of perturbative QCD calculation of 

' hadronic electromaffuetic form factors. 

OO 

cn 

^ ; 1 INTRODUCTION 

0^ ' 

I The applicability of perturbative QCD to exclusive processes at large mo- 

' menta is an interesting research problem. The Brodsky-LePage [||] pQCD 

based factorization has been only partially successful. In this case the 
Qh! process is factorized into a perturbatively calculable hard scattering piece 

and the soft distribution amplitude. The pion electromagnetic form factor 
[^, § at momentum transfer = — Q^, for example, can be written as 



X; F^iQ^) = dxidx2^{x2,Q)H{xi,X2,Q)(j)ixi,Q) (1) 

where (/>(x, Q) are the distribution amplitudes which can be expressed in 
terms of the pion wave function ip{x,kT) as 

= p (fkT^ixM). (2) 

Here x is the longitudinal momentum fraction and kx the transverse momen- 
tum carried by the quark. The factorization is possible provided the external 
photon momentum is much larger than the intrinsic quark transverse mo- 
mentum in which case the kx dependence of the hard scattering H can 
be neglected. 
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The formalism predicts that at large momenta the cross section for ex- 
clusive processes da/dt, where t is the momentum transfer squared, scales 
like up to logs, where n is the total number of elementary partons 

participating in the process. The underlying reason for the power law is 
scale invariance of the fundamental theory. The extra logarithmic depen- 
dence is given by QCD scaling violations. The dominant contribution to 
this scattering arises from the valence quark, since every additional parton 
leads to an additional suppression factor of 1/t. Physically the scattering 
probes the short distance part of the hadron wave function. Dominance by 
the short distance wave functions leads to several predictions such as helicity 
conservation, color transparency [§, ^ etc. 

The successes and failures of this scheme are well known. The predicted 
momentum dependence of exclusive processes, in particular the hadronic 
electromagnetic form factors, have generally been found to be in good agree- 
ment with data. However more detailed dynamical predictions such as he- 
licity conservation in hadron-hadron collisions fail to agree. Calculation of 
electromagnetic form factors using this factorization scheme has been crit- 
icised by several authors. The basic problem is that the momentum scales 
of the exchanged gluons tend to become rather small, and the applicabil- 
ity of pQCD becomes doubtful. The normalization of form factors is largely 
unknown; use of asymptotic distribution amplitudes tends to give small nor- 
malizations compared to data. Form factor magnitudes can be enhanced by 
use of model distribution amplitudes which peak closer to the endpoints, 
namely x ^ 0, 1, which then exacerbates the problem of small internal mo- 
mentum transfers. 

2 THE SUDAKOV FORM FACTOR 

In order to investigate this problem in more detail, Botts and Sterman [^] 
and Li and Sterman developed an alternate factorization which does 
not neglect the kx dependence of the hard scattering. This formalism also 
includes use of a Sudakov form factor. For the case of pion form factor 
the starting point is, 

FniQ"^) = J dXidX2dkTldkT2lp*{x2, kT2,P2)H{xi,X2, Q^, kTi,kT2)ipixi,kTi, Pi), 

(3) 

where it is again assumed that the process is factorizable into hard scattering 
and soft hadronic wave functions ip{x,k'T, P)- The calculation is simplified 
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by dropping the /ct dependence in the quark propagators in hard scattering 
kernel H, in which case only the combination fcyi + ^T2 appears in the 
calculation. The authors Q work in configuration space where this can be 
written as 

f (Pb 

FniQ'^) = J dxidx2-^^-^'P{x2,b,P2, n)H{xi,X2,Q'^,b, fj,)'P{xi,b,Pi,fi), 

where T'{x, b, P, /i) and H{xi, X2, Q^, b, /j,) are the Fourier transforms of the 
wave function and hard scattering respectively; b is conjugate to kxi + kT2, 
fj, is the renormalization scale and Pi, P2 are the initial and final momenta 
of the pion. 

Sudakov form factors are obtained by summing the leading and next to 
leading logarithms using renormalization group (RG) techniques. The wave 
function at small b with a transverse momentum kx cutoff equal to uj = 1/b 
can be approximated by the distribution amplitude (/>(x, 1/6). Large kx 
corrections can be evaluated perturbatively, which result in the Sudakov 
form factor. The final result is given by: 



V{x,b, P, fi) = exp 



-s{x,uj,Q) - s{l - x,LO,Q) - 2 I — 7g(as(/x)) 



X cP{x,l/b) + OiasiLo)). (5) 

where 7^(03 ) is the quark anomalous dimension. The explicit formula for the 
function s{x, u>, fi) is given in Li and Sterman Q. Here lo = 1/b plays the role 
of the factorization scale, above which QCD corrections give the perturbative 
evolution of the wave function P, and below which QCD corrections are 
absorbed into the nonperturbative distribution amplitude (j). For the case of 
the pion, 1/b is the natural choice for this scale. However as discussed below, 
for the proton the relevant scale is not obvious and several possibilities exist 
in the literature. 

The final formula for the form factor, after incorporating the renormal- 
ization group evolution of the hard scattering from the renormalization scale 
/i to t, t = max(y^xiX2Q, 1/b), is given by[^, 

rl roo 

F-i^iQ^) = IOttCf / dxidx2(t){xi)4){x2) / bdbas{t)KQ{./x^Qb) 
Jo JO 

X exp[-5(xi,X2,w,Q)], (6) 
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Figure 1: The Sudakov form factor exp(— S") with = 4 GeV^. For this 
calculation the QCD scale parameter A was taken to 6e 0.1 GeV. 



where 

S{xi,X2, 6, Q) = ^ [s{xi, b, Q) + s(l - Xj, 6, Q)] - 4 —jq{as{fl)). 

The function e~^ is plotted in fig. ||. It cuts off large b regions of the integral 
and hence the calculation is infrared finite, without needing any arbitrary 
infrared cutoff such as a gluon mass. At small b the function has been set 
equal to one. 

The resulting form factor using asymptotic as well as CZ distribu- 
tion amplitudes is shown in fig. |2|. A remarkable fact is that the correct 
asymptotic behavior is seen beyond the scale of about Q = 1 GeV, irre- 
spective of the choice of wave function. In contrast to the Brodsky-LePage 
factorization, the dependence of the hard scattering is not neglected, and 
hence this dependence does not follow trivially. It is rather a detailed 
dynamical prediction of the theory and depends on the relative size of in- 
trinsic kj^ and xiX2Q'^- The prediction is robust, since it is independent of 
the details of the distribution amplitude used. This simple yet important 
point justifies the basic idea of Brodsky-LePage factorization, namely that 
kx can be treated as negligible in the hard scattering. 

We note that the normalization of the theoretical result falls below the 
experimental data for both choices of distribution amplitude. However, the 
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Figure 2: The pion form factor Ft^{Q'^) using the asymptotic (dotted line) 
and the CZ (solid line) distribution amplitudes. The experimental data, 
taken from Ref. 9 with errorbars are also shown. 



large difference between theory and experiment at high momenta should be 



interpreted with caution, since as emphasized by Sterman and Stoler |1C], 
there may be large systematic errors in the experimental extraction of the 
form factor which are not shown in the figure. Further theoretical issues in 
this extraction have been raised in Ref. 

In any event, the theoretical normalization of the form factor is compara- 
tively murky, because it depends on the details of the distribution amplitude. 
Furthermore, the leading order pQCD amplitudes that are practical to cal- 
culate may not give a very reliable estimate of the normalization. One can 
investigate this further by considering the transverse separation cutoff (be) 
dependence of the form factor. This can give an idea about the integration 
regions important for the calculation. We show the be dependence in fig. 
^ as originally discussed in Ref. 0. Based on this plot Li and Sterman 
argue that roughly 50% of the contribution can regarded as perturbative, 
since it is obtained from the region where ag/ir < 0.7. The observation 
also implies that higher order contributions in as are not negligible, and the 
leading order predictions for the normalization of the form factor cannot be 



regarded as accurate. The next to leading order calculation [12] of the pion 
form factor also leads to the same conclusion. 

We are left with the following interesting situation: Although the ba- 
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Figure 3: Dependence of Q'^Ft^{Q'^) on the transverse distance cutoff be with 
the asymptotic (dotted line), and CZ (solid line) distribution amplitudes for 
Q"^ = A GeV^ . The QCD scale parameter A was chosen to be 0.1 GeV for 
this calculation. 



sic Brodsky-LePage factorization is correct, one may need to go to higher 
orders in in order to obtain an accurate prediction for the form factor 
normahzations. However the predicted dependence appears to be quite 
robust, and independent of the theoretical uncertainties such as the choice 
of distribution amphtude. 



3 THE PROTON 

The improved factorization has also been applied to the proton Dirac form 
factor (Q-^)[p!^]. The calculation is considerably more complicated com- 
pared to the pion. Here also it is necessary to use distribution amplitudes 
which peak close to the end points in order to obtain the experimental nor- 
malization of the form factor. In contrast to pion, there is no natural choice 
for the infrared cutoff uj in the Sudakov exponent, due to the presence of 
three quarks and resulting three distances. 

The Sudakov resummation of large logarithms in V leads to 



V{xi,hi,P,fj,) = exp 



s{xi,w, Q)-3 f {osin)) 

1 1 Jw M 



L 1=1 



X (t)[Xi,W) 



(7) 
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where the quark anomalous dimension 7^(0^) = —as/ir in axial gauge gov- 
erns the RG evolution of V. The function (j) is the standard proton distri- 
bution amplitude. The exponent s is given in Ref. 

In equation |^ we use the same infrared parameter uj for all the three 
s{xi,u;,Q) for / = 1,2,3 as well as in the integrals over the anomalous di- 
mension. Earlier Li [^] chose to use different infrared cutoffs bi for each 
exponent s and for each integral involving 7^. As pointed out in this 
choice does not does not suppress the soft divergences from 6/ — > 1/A com- 
pletely, where A is the QCD scale parameter. For example, the divergences 
from 61 1/A, which appear in (j){xi,w) as w — > A, survive as xi — > 0, since 
s{xi,bi,Q) vanishes and s{x2,b2,Q) and s(x3, 63,(5) remain finite. On the 
other hand, w should play the role of the factorization scale, above which 
QCD corrections give the perturbative evolution of the wave function V in 
Eq. (|^, and below which QCD corrections are absorbed into the initial con- 
dition (j). It is then not reasonable to choose the cutoffs bi for the Sudakov 
resummation different from w. 

A modified choice of the cutoffs, w = 1/bmax, bmax = max(6/), I = 1,2, 3, 
was proposed in Ref. |14]. This choice was found to suppress the soft 
enhancements, and the form factor was found to saturate as be 1/A. 
The authors also included a model non-perturbative soft wave function in 
the calculation. Unfortunately, it turned out that the normalization of the 
resulting Q^-Fi was found to be less than half of that of the data for all the 
distribution amplitudes explored |14]. Bolz et al [14| then concluded that 
pQCD is unable to fit the experimental form factor. 

Kundu et al |jl^ reexamined the situation. The group argued that the 
form factor normalization is sensitive to the precise choice of the infrared 
cutoff w. They used w = c/bmax, bmax = max(6/), / = 1,2,3, as the in- 
frared cutoff in the Sudakov exponent, instead of the one used by Bolz et al 
[ 14 1 . The introduction of parameter c is natural from the viewpoint of the 
resummation, since the scale cw, with c of order unity, is as appropriate as 
w Kundu et al |15| find that the calculation is in good agreement with 
data using the King-Sachrajda (KS) [|l6| distribution amplitude and setting 
c= 1.14. 

The choice c = 1.14 can also be motivated physically by considering the 
proton as a quark-diquark type configuration. The diquark constituents are 
the two quarks closest to each other in the transverse plane. Let tityp be 
the distance between the center of mass of the diquark and the remaining 
third quark. Then the infrared cutoff scale u) can be taken to be l/dtyp- We 
choose c such that for a large number of randomly chosen triangles, we get 
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Figure 4: Dependence of Q^Ff on using the KS distribution amplitude 
(c = 1.14, solid line; c = 1, dense-dot line) and the CZ distribution ampli- 
tude (c = 1.14, dashed line; c= 1, dotted line). The experimental data with 
error bars, taken from Ref. 1 7, are also shown. 



for the average (dtyp/bmax) = 1/c. Defining c in such a way, gives c « 1.14. 

The results of the calculation using KS |jT^ and CZ distribution am- 
plitudes and = 1 and 1.14 are shown in fig. It is found that with c = 1.14 
the KS distribution amplitude gives good agreement with data. The be de- 
pendence of the form factor is shown in fig. ^ which shows saturation at 
about be = 0.8/A. The result after including a model nonperturbative soft 
wave function are displayed in fig. ^. Again we find that choosing c of order 
unity gives pQCD calculations in agreement with data. For all choices of the 
distribution amplitude and parameter c, independent of whether the model 
soft wave function is included or not, the dependence of the form factor 
is in good agreement with data. An analogous situation was found for pion 
form factor. 

The natural agreement of dependence of the pQCD calculations is 
in contrast to data fits obtained using soft model |l^, 19 1. In such mod- 
els the dependence depends on the details of the model wave function. 
Soft model predictions at high momentum have a tendency to fall more 
strongly than experimental data. As also noted for the pion, the approxi- 
mate power law behaviour in is not directly implied by the factorization 
and is a detailed dynamical prediction of the calculation. This could have 
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Figure 5: Dependence of Q'^Ff on the cutoff be with the KS distribution 
amplitude for Q"^ = 12 GeV^ (dotted line), Q"^ = 16 GeV^ (dashed line), 
Q2 = 25 GeV^ (dense-dot line), and = 36 GeV^ (solid line). The QCD 
scale parameter A was taken to be 0.2 GeV for this calculation. 

been achieved only if the intrinsic kx were neghgible in the hard scattering 
kernel. Thus the observation of power-law dependence in the data lends 
considerable support to the basic factorization of Brodsky-LePage. Never- 
theless, the relatively small magnitude of internal momentum scales and the 
sensitive dependence of final result on the infrared scale w indicates that 
calculation of normalizations using leading order diagrams is not reliable. 
While higher order contributions are propably non-negligible, there is ev- 
ery reason to believe that the power-law dependence of the calculations is 
robust. 

We have reviewed the current status of pQCD calculations of hadronic 
electromagnetic form factors. We argue that the normalization of the form 
factor cannot be predicted reliably by a leading order calculation in Ug- De- 
tailed calculations including the soft kx dependence, however, support the 
basic factorization scheme. One finds the correct asymptotic evolution 
of the form factors for Q as small as 2-3 GeV, independent of the choice of 
distribution amplitude and other theoretical uncertainties. Hence we con- 
clude that agreement of quark counting scaling predictions is not accidental 
but is well supported by detailed dynamical calculations. 

PJ and BK thank the staff of ICTP, Trieste, for hospitality during a 
visit where this paper was written. This work was supported by BRNS 
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Figure 6: dependence of Q'^Ff including a soft model wave function with 
< kj^ >= 0.271 GeV^. The four upper curves at = 35 GeV^ use the 
KS distribution amplitude with the infrared parameter c = 1.4 (solid curve), 
c = 1.2 (dashed curve) and c = 1 (dashed- dotted curve). The dotted curve 
shows the result without including model soft wave function with c = 1. The 
two lower curves at = 35 GeV^ use CZ distribution amplitude with c = 1 . 
The lowest curve includes the soft model wave function, whereas the upper 
curve does not. The experimental data with error bars [17] are also shown. 
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